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. . . We show that the long wavelength modes of a field become classical during a second order phase 

• transition because of the interaction with the short wavelength modes of the field. In a massive scalar 

' frield model the number and thermal states of long wavelength modes, whose Wigner functions are 

^— , sharply peaked around the classical trajectories during the phase transition, exhibit only classical 

correlation without achieving quantum decoherence. In a linearly coupled scalar field model, the 
^ ' long wavelength modes are shown to efi'ectively achieve quantum decoherence because of the mode- 

Q , mixing. Finally we define a quantal ordering parameter that is linear in the field variable and 

■ satisfies the classical field equation. 

CS| I PACS numbers: 03.70.+k, ll.lO.Wx, 05.70.Fh, 11.30.Qc 

I ! 
(N ■ 

> : 

, I. INTRODUCTION 

(N ■ 

A system becomes classical when it recovers classical correlation and loses quantum coherence, i.e. decoheres. The 
\ present Universe, for instance, is believed to experience a quantum-to-classical transition at a certain stage of its 
. evolution from its early quantum state. One of such decoherence mechanisms is the interaction of the system with an 

■ environment ^, Q| . The environment-induced decoherence mechanism explains correctly the density perturbation 
, necessary for structure formation of the present Universe |^ . A phase transition is another interesting phenomenon 

O i' in which one observes the quantum-to-classical transition. For instance, the system of the long wavelength modes of 
I ] a self-interacting scalar field become classical through the interaction with the environment of the short wavelength 
modes |6). 

(-H I Another aspects of classicality is classical correlation. The parametric interaction of an open system provides 
• • . classical correlation without any direct interaction with the environment. The coupling constants (parameters) of the 
open system change explicitly in time. A quantum field model in which the mass parameter changes sign during a 
quenched second order phase transition provides such an open system. In this model the dynamical evolution of phase 
%^ ' transition is mainly described by a classical order parameter, whose quantum fluctuations are necessary for domains 
or topological defects 0. In a previous paper |^ (hereafter referred to I), we introduced a quantum phase transition 
model without any classical order parameter in which a certain symmetric quantum state drives the phase transition. 
A question is then raised how the quantum phase transition exhibits classical features. 

In the field model a classical background field is introduced as the order parameter, around which quantum field 
fluctuates [Q. The expectation value of the field, however, vanishes globally from the symmetry (parity) argument 
[^. There is still a room for explanation of the order parameter in quantum theory: the condensation of the field 
into a coherent state having a nonzero expectation value. However, in Ref. I this initial quantum state is shown 
to be either a coherent or a coherent-thermal state. If the system starts initially from an exactly symmetric state 
such as the Gaussian vacuum or thermal equilibrium, the quantum law yields zero expectation value throughout the 
phase transition. On the other hand, in Ref. I another possibility is discussed that the phase transition may proceed 
quantum mechanically in the symmetric quantum state, such as the Gaussian vacuum, number, and thermal states, 
even without globally forming a coherent condensate. 

In this paper we shall show how classicality emerges from the symmetric quantum evolution of phase transition. Our 
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model is a massive scalar field, the mass of which changes sign to mimick a quenched second order phase transition. 
This model initially has been used to study the behavior of an inflaton during the slow-rollover in the new inflationary 
scenario [|j and also has been employed to describe quantum processes of the phase transition during the spinodal 
instability regime in Ref. I. It is known that the Gaussian state of the model obeys a classical probability distribution 
. In this paper we shall apply the quantitative measure of both classical correlation and quantum decoherence in 
Ref. []lo| to the quantum phase transition. For that purpose we derive the density matrices and Wigncr functions for 
the Gaussian vacuum, number state, and thermal equilibrium. 

It is found that the density matrix and Wigner function manifestly exhibit classical correlation for the long wave- 
length (soft) modes that grow exponentially, whereas both long and short wavelength (hard) modes maintain the 
same initial quantum coherence. Hence the massive scalar field does not achieve classicality in the genuine sense. 
To show how quantum decoherence occurs during the phase transition, we consider an analytically solvable model, 
motivated by a self-interacting scalar field, in which a long wavelength mode is linearly coupled to a short wavelength 
mode. It is shown that the mode-mixing (coupling) between the long and short wavelength modes drastically reduces 
the degree of quantum coherence of the long wavelength mode during the phase transition. Thus the long wavelength 
mode becomes completely classical from the view point of both classical correlation and quantum decoherence. It is 
further suggested that a quantal quantity can be defined out of long wavelength modes to play a similar role of the 
order parameter. 

The organization of this paper is as follows. In Sec. II, we derive the density matrices and the Wigncr functions 
of a massive scalar field in the Gaussian vacuum, number state, and thermal equilibrium. In Sec. Ill, using the 
quantitative measure of classical correlation and quantum decoherence, we show that the unstable long wavelength 
modes exhibit classical correlation, whereas the short wavelength modes maintain their initial quantum coherence. In 
Sec. IV, we show how a long wavelength mode can achieve quantum decoherence during the phase transition through 
a linear coupling with a short wavelength mode. In Sec. V, we define a quantal order parameter that is linear in the 
field variable and satisfies the classical field equation. 



II. DENSITY MATRIX AND WIGNER FUNCTION 



As the first model for the second order phase transition, we consider a real massive scalar field with the Lagrangian 
density 



/: = i[02-(v</))2]-im2(^)0^ 



(1) 



where the mass is assumed to change sign during the phase transition. Following Refs. I and the action (|^), 
upon a suitable mode-decomposition, leads to the Hamiltonian 



m = E 



where 



and a denotes the Fourier mode defined by 



^i^\t) = -[Mt) + <t>Mt)], A~\t) = '-[Mt) - 0-k(i)]- 



(2) 



(3) 



(4) 



The Fourier modes (Q) are Hermitian because <pl. = (p-^. Thus the Hamiltonian (||) consists of the infinite number of 
decoupled, time-dependent oscillators. In quantum field approach the fundamental law is the functional Schrodinger 
equation (in units of h = k = 1) 



i^vE'(0,<)-i?(t)*(0,t). 



(5) 



As every mode is decoupled from each other, the wave functional to Eq. (^ is now given by the product 

vl/(0,t)=[]vl;„(0„,i) (6) 
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of the wave functions of each Schrodinger equation 

d 



at 



(7) 



The Fock space consisting of exact quantum states for a time-dependent oscillator was first constructed by Lewis and 
Riesenfeld fT^ (for the references on many different methods, see [p^). 

Now the quantum evolution of the Hamiltonian is reduced to that of individual oscillators in Eq. (^). To find 
each Fock space, we follow the so-called Liouville-von Neumann (LvN) approach developed in Refs. I and pi] , |li 
the LvN approach one finds a pair of operators for each mode of the Hamiltonian (0) : 



that satisfy the quantum LvN equation 



d r ait) 



a{t) 
at(t) 



Equation 



'dt \ a^t) 

leads to the classical equation of motion for a complex ifa- 

(pa{t) + Ujl{t)ipo,{t) = 0. 

By requiring the Wronskian condition 

~ <y5a(iX(i) = j, 

one can make haify and a]j,(t) satisfy the standard commutation relation at any time 



In 

(8) 

(9) 

(10) 
(11) 
(12) 



In the LvN approach the exact quantum state of the time-dependent oscillator is determined, up to some time- 
dependent factor, by the eigenstate of a Hermitian operator satisfying Eq. (|^). For instance, the vacuum state is 
the zero-particle state aQ.(i)|0,t) = according to the standard quantum mechanics. Hence, as the phase transition 
proceeds, the time-dependent vacuum state at a late time is a squeezed state of the initial vacuum state [ p^ p^ . 
As both operators in Eq. (^) already satisfy Eq. (||), we choose each number operator as 



7V„(t) =4(t)a„(0, 
and construct the Fock space consisting of the time-dependent number state 

Na(t)\na,t) = na\na,t). 



Then the vacuum state is given by 
and the number state by 



aa(i)|Oa,i> =0, 



In Rcf. I the number state has the coordinate representation 

/ \ 1/4 ^ / X 

^nj^i-a,*) 



v/2 



exp 



2 



(13) 
(14) 
(15) 

(16) 
(17) 



where i?„^ is a Hermite polynomial. From the definition p^i,{x' ^x) = {x'\^){^!\x) , we find the density matrix for the 
number state (p7[): 



1 



1/2 



X exp 



1 



+ ln('^* '^a)(?!)a,C0Q,A 



(18) 
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where 



Once given the density matrix (|l^), one easily finds the Wigner function jl^ 

f /■+°° 

Pa(0a,7ra) = - / - ^ |/5a W 1 0a + 2/) 6^""°^ . 

I" J-oc 

Substituting Eq. ( p^ ) into Eq. (pO|), we obtain the Wigner function for Ha 

P«„(0a,7r„) = -(-l)"°L„„(8(^>„#„)e-4'^>°^°, 

TT 

where is a Laguerre polynomial and 



TTa - < ^in(v3;!;.^a)^/^ 



H 2 



1 



8(^:;^a)2'""- 



(19) 



(20) 



(21) 



(22) 



We now turn to the density matrix and Wigner function for the thermal state. In Refs. I and [ pj] |, the following 
density operator is used for each mode 



PTa{t) — 2 sinh 



f3uj. 



(23) 



because Na{t) already satisfies the quantum LvN equation. Here the free parameters (3 and UaA are to be identified 
with the temperature and frequency, respectively, of the initial thermal equilibrium. From the density matrix of Ref. 
I 



PTa(0a,0a,O 



tanh(/?Wc[,i/2) 



1/2 



exp 



X exp 



1 



2Va'^a 



tanh 



we obtain the Wigner function 



1 



fTa(0a,7i'a) = - tanh 



[ 2 



exp 



-4 tanh 



(24) 



(25) 



Also, the density matrix ( |18[ ) and Wigner function ( pi p for the Gaussian vacuum state (n^ = 0) are obtained by 
taking the zero temperature limit /3 — > oo of Eqs. (^J) an d (|25| ). 

The nature of the number state ( |l7| ) and thermal state (^3|) can also be understood from the Wigner functions ( [2l| ) 
and (|2^). In particular, the 1-a contour of the Wigner function contains some useful kinematical information of the 
quantum state in phase space. The 1-a contour of the thermal state is defined by 



tanh 



V 2 



>A(p^(paHa = 1. 



(26) 



Here and hereafter the 1-a contour of the Gaussian vacuum state is obtained by taking limit tanh(/?a-'Q,^i/2) = 1. As 
Eq. ( p6| ) describes a rotated ellipse in the phase space {TTa,<j>a), we can find the major and minor axes through the 
transformation 



Tin = COS Un TTry — SlU ( 



= sinOaTTa + COS ( 



(27) 
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where 



tan 29 n 



2Ca 

Ba — Ba 





tanh 


/ o \ 1 
\ 2 )\ 




a = < 










tanh 


( (iiOaA 




1 


K 2 )\ 







Then Eq. (Eq) is written in a canonical form 



(28) 



where 



1 ^ Bt+B- 
A± 1^ 2 

Physically, the constant area of the ellipse given by 

TT 



± 



B+-B. 



vAJAo 



coth 



a a 



cos( 



(29) 



(30) 



(31) 



implies that the phase transition only squeezes the initial quantum state of the massive scalar field unless there is a 
mode- mixing (coupling). The entropy of each mode given by 



Sa ~ In Aa = In coth 



\ 2 r 



(32) 



is also constant, i.e. isentropic, and vanishes for the vacuum state as expected. Without coarse graining the entropy 
of the massive scalar field does not change. 



III. CLASSICAL CORRELATION OF LONG WAVELENGTH MODES 



We now study classical correlation and quantum decoherence of the density matrix ([T^) or ( p^ and the Wigner 
function (21) or (p5|). The Wigner function ( pl| ) oscillates for an excited state and is not always positive definite, so it 
does not describe a true probability distribution in phase space. On the other hand, the Wigner function is positive 
definite for the Gaussian vacuum and thermal state. The thermal state is particularly of physical interest in the phase 
transition because it is mostly likely that the system starts from a thermal equilibrium. 

A quantum system achieves quantum decoherence when the interference between classical trajectories is lost, and 
it recovers classical correlation when the Wigner function is peaked along a classical trajectory. The coherence length 
Ix of the density matrix ([l^) or (^4|) written in the form 



tanh 



PTa, 



(/3c^„,z/2)}' 



1/2 



27r<^* V3q 



exp 



Tee, A' 



M(Pa,C(Pa,A 



(33) 



is roughly determined by the width of 0^ ^ : 




1/2 



(34) 
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Here the vacuum state result is ta.nh{(3uja,i/'2) = 1- Hence the coherence length increases or decreases depending on 
whether ipa grows or decays. A large coherence length implies a large degree of quantum interference, so quantum 
decoherence is conditioned by the small magnitude of Ix- 

More rigorously, the representation-independent measure of classical correlation of Ref. jl^ is given by 



Sec = 



^ a,C^ a, A 
1 



(35) 



Similarly, the measure of quantum decoherence , which is determined by the ratio of the width of the off-diagonal 
element to that of the diagonal element, is given by 



3qd 




(36) 



This means that the massive scalar field preserves the same initial quantum coherence throughout the phase transition. 

The quantum system recovers classical correlation when 6cc 1 and loses quantum coherence when Jqd ^ 1. So 
the system becomes classically correlated only when ipa grows during the evolution. In the opposite case of decaying 
ifa, the degree of classical correlation increases compared to that of quantum coherence, but the system retains 
quantum coherence. To see how the classical trajectory is recovered when dec ^ 1, we look into the Wigner function 
(2^). Each contour of the Wigner function, which is quadratic ttq, and (j)a, depicts an ellipse in phase space. Suppose 
that both the classical field 4>a and auxiliary field ipa grow exponentially during the phase transition: 



fa (t) 



then the classical trajectory obeys 



(37) 



(38) 



On the other hand, the last term in Eq. (|2^) is exponentially suppressed compared to the first two terms and the 
overall coefficient of the exponent increases exponentially, so the Wigner function ( p5| ) is sharply peaked around a 
trajectory in phase space: 



0Q — fa{t)<j)a- 



(39) 



Indeed the Wigner function is sharply peaked around the classical trajectory. 

Now we apply the criterion on classicality to the second order phase transition via an instantaneous quench, in 
which the mass changes sign as 



ml t < 0, 
i-my, i > 0. 

Before the quench, the most general solution to Eq. ( |l0[ ) is given by 

1 



(40) 



fo 



cosh r^e 



-iuJaAt 



sinh r„ e 



iUa^it 



where 



k2. 



(41) 



(42) 
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The system is stable and has the minimum expectation value when Tq, = 0, corresponding to the true vacuum state. 
As rQ,(^ 0) represents a squeezing parameter of the true vacuum state, the vacuum state in Eq. ( |l6| ) is the one- 
parameter squeezed vacua After the quench, the solution for each unstable long wavelength mode matches the 
initial solution (^) continuously at the onset of the quench {t — 0) and is given by 



fa,Uf 



(coshrQ,-|-e sinhra) cosh{ujajt) 



-i(coshrQ, 



sinhra) _ sinh(a3Q, ft) 



where 



mj — 



(43) 



(44) 



Whereas the solution for each stable short wavelength mode is given by 



'Pa,Sf 



(cosh 



sinh Va) cos{uja, ft) 



-i{coshra — e * ° sinhro.) sin(ci;a ft) 

UJaJ 



where 



k-^ — rrij. 



Thus, after the quench, the amplitude square of the long wavelength mode increases exponentially as 

1 

■J 



fa UfVa.Uf = TTz [cosh(2rQ) cosh(2(DQ./t) + sinJa sinh(2rQ) sinh(2a)Q./i)]. 



(45) 



(46) 



(47) 



Therefore, the measure of classical correlation (35) for the long wavelength modes decreases exponentially to zero as 
the spinodal instability continues. Also, the Wigner function (^ij) is sharply peaked around the classical trajectory 
TTa — Cjaj4>a- On the Contrary, the short wavelength modes have oscillating, bounded, amplitudes 



[cosh(2rQ) cos{2uJaj-t) + sin^Q sinh(2rQ) sm{2uj a, ft)]. 



(48) 



Thus the measure of classical correlation ( p5| ) cannot become small for the short wavelength modes. But the measure 
of quantum decoherence ( ^6|) of the long and short wavelength modes remains constant. 

In summary, the unstable long wavelength modes of the massive scalar field recover classical correlation during the 
phase transition whereas the stable short wavelength modes do not. However, both the long and short wavelength 
modes retain the same quantum coherence throughout the phase transition. 



IV. QUANTUM DECOHERENCE DUE TO MODE-MIXING 



The long wavelength modes of the massive scalar field model in Sec. Ill gain classical correlation during the 
second order phase transition without achieving quantum decoherence. A more realistic model of the quantum phase 
transition is provided by $^ theory with the Lagrangian density 



/: = i[<^^-(v<^)2]-im^(^)</>^-A^^ 



(49) 



where m?{t) is assumed to change sign as in Eq. 
Hamiltonian obtained from Eq. (^9|) takes the form~~||I 



Upon decomposing the field into the Fourier mode (^, the 

(50) 



4! 
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where lo\ is given by Eq. (^). The last coupled anharmonic terms prohibit any analytical solution. In the Gaussian 
approximation of the Hartree-Fock or Liouville-von Neumann method ||] , the number state of each mode is still 
given by Eq. (|T7|), where (/Sq now obeys the mean-field equation 



(pa (t) 



=0. 



(51) 



In the Gaussian approximation the term c/)'^ affects only the auxiliary variable ipa through Eq. (pl[), but does not 
change the form of the state (17). Therefore, the Gaussian approximation for the self-interacting scalar field (^ ) does 
not lead to quantum decoherence. As the short wavelength modes provide an environment to the long wavelength 
modes, the mode-mixing is expected to lead to quantum decoherence. 

To analytically study the mode-mixing effect between a long and short wavelength mode, we consider an exactly 
solvable model 



(52) 



Here the subscript S and H denotes the long and short wavelength mode, respectively. As Lojj and tojj depend on time 
through m?{t), the orthogonal transformation leading to normal modes depends time explicitly. Hence the product 
of the instantaneous eigenstates of each normal mode is not an exact quantum state of the original Hamiltonian (|5^) . 
The Schrodinger equation has a Gaussian state of the form 



*o( 



,0ff,t) = 7V(t)exp -{As{t)(l)l + \B{t)c^s(l)H + AH{t)4>]j} 



(53) 



where N is the normalization constant, which is not of much concern to this paper. The coefficients determined by 
the Schrodinger equation are 



where 



u{t) 



As{t) = - 
B{t) = i 



'it) 



2u*{ty 
2v*{ty 

Ju*{t)v*{t) 




0, 



We solve Eqs. ( p5| ) and ( |56| ) separately for two cases: before and after the phase transition. 

First, before the phase transition, the frequencies take the constant values ojgit) = uig ^ and uj'jf{t) — 



solutions to Eqs. (p5[) and (pq) are found to be 



'HA- 



(54) 

(55) 
(56) 

The 



v{t) 



1 



1 



(57) 



where 
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Hence the Gaussian state is given by 



and the reduced density matrix for the long wavelength mode by 



Pr('/'s>s) = 



+ 00 



N*N, I — exp 



Finally, from the coefficients of the density matrix 

A2 



s.c 



the measure of quantum decoherence is given by 



rs,A = ^s, ^s,M = 0, 



(59) 



(60) 



Sqd = - 



1 Ts. 



c 



2 V rs,A 



and that of classical correlation by 



'Jcc 



p2 p2 



(61) 



(62) 



Therefore, quantum decoherence is achieved by the mode- mixing, but the characteristic behavior of classical correlation 
does not change. 

Second, after the phase transition, Eqs. (pq) and (p6h become 



Hit) ^ 

m 



u{t)v{t) 

Sn{t)v{t) \ 
u{t)v{t) I 



2n 



u{t) = 
v{t) = 0. 



In the weak coupling limit M <^ Cjsj,liJhJi the approximate solutions are found to be 



(63) 



(64) 



where 





1 ^sj 4 























A^ / 




- i^Hj 










( ^sj - 


h iujHj 









(65) 



These solutions are accurate at late times provided that Re(ri5)i 3> 1. Then, the reduced density matrix is given by 
Pr(0s, 05) = ^*^y'^^ ^^^ exp[-rg,c0s,c - rs,A0|,A - ^s.m(I)s,c(I>s.a] , (66) 
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where 



S,A 



= 2 



S,M 



ReAs 



ReAc 



ARoAh 

AReAn 

A2(ReB)(ImB) 
2ReAH 



To order A^, the coefficients are approximated by 



s,c 



A2 



1 



X^UJHJ 



r 



S,M 



sj 



,2 



2u) 



SJ 



16 Sujjj 



2^lf 



'UJHJ 



(67) 
(68) 
(69) 

(70) 
(71) 
(72) 



Therefore, to A^ order, the measure of quantum decoherence is given by 



2\j 2Qs,fi^lj+u;j,jf 



and that of classical correlation by 



5A6 



3CC 



\UJS.J - UJH.jl 



4 C^lfiuj J + ajj^jY 



(73) 



(74) 



A few comments are in order. First, quantum decoherence is observed even for the quadratic Hamiltonian (|5£ 
It is a consequence of the mode-mixing and does not exclusively pertain to the nonlinearity of the system. For the 
quadratic system with constant parameters, the parameters introduced in Ref. B as 




tan(26') 



A 



9 9 ' 



where 



.|,)2+A2 



(75) 



(76) 



yield another expression for quantum decoherence 



2 y cosh — sinh ricos^{29) 



(77) 



Obviously, there is no quantum decoherence, Sqd — 1/2, for the zero mode-mixing X = 9 = 0. In the limiting case of 
weak coupling A <C — a;||, the mixing angle is small {9 « 0) and the measure of quantum decoherence becomes 
Sqd ~ 1/2, which is consistent with the result (|6l|). Whereas, for the two identical oscillators us — uh with the 
maximum mixing angle a = 7r/4, the measure is simply given by 



Sqd = 



1 

2 cosh rj 



(78) 
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where 

/ 2a)| + A 

4 



Second, it is shown that the unstable system efficiently achieves quantum decoherence through the coupling with an 
environment. This effect is also observed in the Bateman or the Feshbach-Tikochinsky oscillator. The Bateman or 
Feshbach-Tikochinsky oscillator is a conserved two-oscillator system which consists of a damped and an amplified 
oscillator [0, |l^. The amplified oscillator has an exponentially increasing amplitude at the expense of the damped 
oscillator and reminds us of the unstable mode of the model (p^. It is found that the amplified oscillator shows a 
similar quantum decoherence [^ . Thus it may be inferred that even in quadratic systems quantum decoherence is a 
consequence of the mode-mixing (coupling) of systems to the environment. 

V. QUANTAL ORDER PARAMETER 

In the previous section it is found that the unstable long wavelength modes recover their classical correlation but 
the stable short wavelength modes retain the quantum coherence during the quench. Thus the parametric interaction 
via the quench makes the quantum phase transition exhibit classicality. As some, though not all, of the essential 
features of the second order phase transition can also be explained by classical theory, it would be interesting to find 
a quantal quantity that behaves like a classical order parameter. 

In quantum theory the order parameter is frequently defined by ((/)(i,x)), the mean value of the field. The mean 
value of a coherent state satisfies the classical field equation with quantum correction from fluctuations. However, 
the thermal state or Gaussian vacuum studied in this paper is symmetric, so the expectation value of the fleld modes 
(field) with respect to this state vanishes: 



MLpT,a) = 0. (80) 



The nonzero expectation values come from the quadratic moment of the field or field modes. One candidate is the 
two-point function (0(t, x)</)(t, y)). It satisfies the classical field equation in the case of a massive scalar field but not 
in the case of a self-interacting scalar field. As the two-point function is quadratic in the field variable, we wish to 
define a quantity that is linear in the field variable and satisfies the classical field equation. 
We note that the quadratic moment of the field mode also has the nonzero expectation value 



12 \ 



T = ¥?*(PaC0th|^^^J. (81) 



The classical equation (|T^) for a complex (fa leads to the equation 



where is the amplitude of complex tpa- 



la (t) + u;l {t)U (t) - - 0, (82) 

^cit) ^ Ut)e''-^'^ ■ (83) 

The last term in Eq. (^2[), which originates from quantization rule and guarantees the Heisenberg uncertainty principle, 
can be neglected for the long wavelength modes because increases exponentially during the quench. Therefore, it 
approximately satisfies 

Ut)+u;l{t)Ut)^0. (84) 
As each ^ke'*' '' satisfies the classical field equation, it follows that the quantal quantity defined by 



<^o(t,x) = / TTT^Ut) 



/>k<7n f 




Jk^O 


(27r)3 


pk<mf 




/fc=0 


(27r)3 



^^)Ttanh(^^ 



-I 1/2 



e*'^-'' (85) 
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also satisfies the classical field equation 

^(i,x)- V20o(i,x)+m2(t)0o(i,x) ~0. (86) 

It should be remarked that (jfo is linear in the field variable and is defined out of the long wavelength modes, which 
become classically correlated during the quench. In this sense 4>o carries many features of the classical order parameter. 

VI. CONCLUSION 

We studied classical correlation and quantum decoherence in quantum phase transition. For the initial state such as 
the Gaussian vacuum or number state or thermal equilibrium, the expectation value of the field vanishes throughout 
the phase transition. Hence it cannot be used as the order parameter. The phase transition in such a symmetric state, 
however, can be explained entirely within the framework of quantum field theory. In the massive scalar field model the 
long wavelength modes become unstable and grow exponentially during the quench. The quantitative measure from 
the density matrix shows classical correlation of the long wavelength modes during the phase transition. Likewise, 
the Wigner functions for these modes are sharply peaked around their classical trajectories and thus confirm classical 
correlation. Whereas, the short wavelength modes are stable throughout the quench process and retain the quantum 
coherence. 

Classicality is, in a strict sense, conditioned not only by the recovery of classical correlation but also by the loss of 
quantum coherence. However the massive scalar model does not achieve completely quantum decoherence, because the 
measure of quantum decoherence keeps the value determined by the initial Gaussian vacuum or thermal equilibrium. 
Hence the system recovers classicality partially, depending on the relative magnitude between the measures of classical 
correlation and quantum decoherence. To achieve quantum decoherence it is necessary for the long wavelength modes 
to couple to an environment. The environment-induced decoherence was studied numerically in the self-interacting 
scalar field model, where the long wavelength modes are coupled to the short wavelength modes |2^. However, it 
is very hard to find any analytical solution of the self-interacting scalar model. The Gaussian approximation of the 
Hartree-Fock Q or the Liouville-von Neumann method Q does not include the mode-mixing effect appropriately to 
explain quantum decoherence. 

In this paper, to analytically study the mode-mixing effect, we turned on a linear coupling between the long and 
short wavelength modes of the massive scalar field. This exactly solvable model shows that the long wavelength modes 
indeed decohere because of the mode-mixing with the short wavelength modes. Therefore, classicality emerges during 
the phase transition from the unstable long wavelength modes coupled to the stable short wavelength modes. These 
long wavelength modes behave classically and provide a quantal quantity that is linear in the field variable, satisfies 
the classical field equation, and thus behaves like a classical order parameter. 

After this paper was completed, we were informed of Ref . pH | , which also shows quantum decoherence of the long 
wavelength modes through the nonlinear coupling to the short wavelength modes in the self-interacting scalar field 
and confirms the numerical result |^. However, any analytical solution in a closed form has not been found for the 
self-interacting scalar field. At best, the analytical solution can be found by perturbatively including the mode-mixing 
terms |p^ , which is beyond the scope of this paper. The non-Gaussian effect from the mode-mixing is expected to 
result in quantum decoherence. As quantum decoherence is largely a consequence of mode- mixing (coupling), our 
exactly solvable model may shed some light in understanding the mechanism of classicality. 
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